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ABSTRACT 



We present an up-to-date analysis of the constraints imposed bythe precision data on the (CP- 
conserving) Two-Higgs-Doublet Model of type II, with emphasis on the possible existence of 
very light neutral (pseudo) scalar Higgs boson with mass below 20-30 GeV. We show that even 
in the presence of such light particles, the 2HDM(II) can describe the electroweak data with 
precision comparable to that given by the SM. Particularly interesting lower limits on the mass 
of the lighter neutral CP— even scalar are obtained in the scenario with a light CP— odd 
Higgs boson and large tan/3. 
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1 Introduction 



The Standard Model (SM) of electroweak interactions is in very good agreement with the elec- 
troweak precision data collected at LEP and SLAC experiments This fact strongly supports 
the idea of the spontaneous breaking of the underlying SUl{2) x UyiX) gauge symmetry. Yet, 
the actual mechanism of symmetry breaking still remains unexplored. The Higgs boson pre- 
dicted by the minimal model of electroweak symmetry breaking, the SM, has not been found 
up to now. Only the lower limit on its mass of ~ 90 GeV is set by the unsuccessful search 
at LEP. While the extreme simplicity of the Higgs sector of the SM is theoretically appealing, 
there exist many of its extensions that lead to different phenomenology (more physical Higgs 
particles) and which also should be tested (or constrained) experimentally. 

The simplest such extension is the well-known Two-Higgs-Doublet model (2HDM). It exists 
in several distinct versions, of which we want to consider in this article the one that is called 
Model II in its CP— conserving version (we briefly recall its structure in the next section) .0 
One interesting question that arises in the context of such an extension of the SM is what 
the available experimental limits are on the masses of the Higgs bosons predicted in such a 
model. In Section 3 we will recall the arguments that, in the framework of the considered 
version of 2HDM, the direct searches do not exclude the existence of very light neutral scalar 
or pseudoscalar Higgs particles. In Section 4 we show that the existence of such light Higgs 
bosons is notexcluded by the electroweak precision data either. An analysis of the impact of the 
precision data on the 2HDM(II) was already performed in the past ^ (for the formalism and 
early investigations, see also |Q) but concentrated mainly on the possibility of improving the 
prediction for Rb = r{Z^ —>■ hh)/T{Z^ hadrons) (which at that time seemed to be required 
by the data) and on improving the 2HDM(II) global fit to the data with respect to the fit 
given by the SM. Since then, the experimental situation has evolved significantly. In particular 
the measurement of no longer shows any statistically significant deviation from the value 
predicted by the SM [0. Also there are changes both in the experimental measurement and 
theoretical computation of the 6 — 57 decay rate, which was crucial in the analysis performed 
in ref. 0. More recently a partial analysis of the constraints imposed on the 2HDM(II) by 
various measurements was also attempted in Here we present an up-to-date analysis of the 
constraints that the precision data impose on the 2HDM(II), with emphasis on the possible 
existence of light neutral (pseudo) scalar Higgs boson. We show that even in the presence 
of such light particles, the 2HDM(II) can describe the electroweak data with the precision 
comparable to that given by the SM. In this case some interesting global limits on the model 
can be obtained. Finally in Section 5 we summarize our results and briefly comment on the 
other ways the existence of light Higgs bosons are or can be constrained by other experimental 
data. 



The Higgs sector of the Minimal Supersymmetric Standard Model has precisely the structure of Model II, 
but with additional constraints imposed on the quartic couplings. However, in the following we will consider 
the general Model II in a regime in which it cannot be regarded as a low-energy approximation of the MSSM 
with heavy sparticles. 
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2 Two-Higgs-Doublet extension of the SM - Model II 



The muhidoublet extensions of the SM are distinguished by their virtue of not introducing 
corrections to the p parameter at tree level. The minimal extension of the SM consists of 
two doublets. The requirement of the absence at tree level of the flavour- changing neutral 
currents puts restrictions on how the two scalar doublets of the general 2HDM can couple to 
fermions. In Model II, one Higgs doublet (denoted by $i) couples only to leptons and down- 
type quarks, whereas the other doublet ($2) couples only to up-type quarks. After spontaneous 
symmetry breaking, both doublets acquire vacuum expectation values Vi and f2, respectively 



— =tan/5. (1) 

Vl 



with V = J vf + V2 fixed by Mw, and 



With two complex Higgs doublets, the 2HDM predicts the existence of five physical scalars: 
neutral h^, and and charged H^. In the more restrictive scenario (which we are going to 
discuss), with CP symmetry conserved by the Higgs potential, and are CP— even mix- 
tures of the neutral components of the doublets (the mixing being parametrized by the angle a) 
whereas is CP— odd. Thus, in the CP— conserving version, the Higgs sector is parametrized 
by four masses M/j, Mh (by definition Mh < M^), Ma and Mu+ and three dimensionless 
parameters tan/3, a and A5 (for definitions see for instance, 0). The first two dimensionless 
parameters are very important for the phenomenology of the Higgs sector as they determine 
the couplings of the physical Higgs bosons to fermions and gauge bosons: the couplings of the 
scalars to the down- and up-type quarks are given by the SM couplings multiplied by the factors 
(see for instance, |]^, ^) 

h^bb : = sm(p — a) — tan P cos(/? — a) 

cos p 

— COS (y, 

h^tt : = sin(/5 — a) -|- cot (3 cos(/3 — a) 

sin p 

_ cos (y 

H bb : = cos(/5 — a) -|- tan /5 sin(/3 — a) 

cos p 

HHt : = cos{/3 -a)- cot /3 sm{(3 - a). (2) 

sin p 

The Feynman rules for the CP— odd scalar couplings to fermions are given by the SM rules for 
/i^^j times the factors: 

A%b: -i-fHan(3, AHt : - cot /5. (3) 

Important for the direct Higgs boson search at LEP, couplings Z^Z^h^ and Z^Z^H^ are given 
by the corresponding Standard Model coupling Z^Z^h%j^.j modified by the factors: 

Z°Z°/i°: sin(/?-«), Z^'Z^'H'': cos(/?-«), (4) 

whereas the couplings of Z^ to A^h^ and A^H^ pairs are instead proportional to 

ZM°/i°: cos(/?-a), Z^A^H^: sm{(3-a). (5) 
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Figure 1: 95% CL. limits on sin^(/5 — a) from the Higgs boson search at LEP as a function of 
Mh. Lines 1 and 2 are the OPAL results obtained assuming the Higgs boson decay branching 
ratios as in the SM and 100% into hadrons, respectively [0. Line 3 is the result of L3 |12|. 



It follows that, in the limit sin(/3 — a) = 1, the lighter CP— even neutral Higgs boson has 
precisely the couplings of the Standard Model Higgs and becomes indistinguishable from it. 
Finally, couplings of the charged Higgs scalar to fermions, e.g. btH~ vertex, are given by 
expressions like: 

^ [mtCot/5(l +75) +mbtan/3(l -75)] . 

It is clear that for tan/3 close to zero, scalar couplings to tt pairs are strongly enhanced with 
respect to to the SM case, whereas for tan/3 large (> 10) the couplings to \k are enhanced. 
Requirement of perturbativity of both types of couplings restricts, therefore, tan /? values to 
the range 0: 

0.3 < tan /3 < 130. (6) 
Outside this range perturbativity is lost and no firm prediction can be obtained from the model. 



3 Constraints on the 2HDM(II) from direct Higgs boson 
search at LEP 



Here we briefly review the constraints on the 2HDM(II) imposed by direct Higgs boson search 
at LEP [0], which will be taken into account in performing the fits to the electroweak data in 



Sec. 4. At the LEP collider, the Higgs boson search is based mainly on the following processes: 



1. the Bjorken process e~^e~ — >■ Z^*h^{H^) 

2. the associated Higgs boson pair production e^e~ A^hP{H^) 

3. the Yukawa process e^e~ bb ^ bbA^{h^) 

4. e+e" ^ H+H-. 

Because of the structure of the Z^Z^hP and Z^A^h^ couplings (|D,(H) the processes e"'"e~ 
Z^*hP and e^e~ —>■ A^h^ are complementary to each other, provided they are simultaneously 
kinematically allowed. In the SM, or in the 2HDM(II) in the limit sin(/5 — a) = 1, the non- 
observation of the Bjorken process at LEP sets a lower limit on Mh^j^j of ~ 90 GeV (at 95% 
C.L.) [1^. In the general case, the same data put an upper limit on the factor sin^(/? — a) 
as a function of the lighter Higgs boson mass (see Fig. M. The combined analysis of the 



complementary channels 1 and 2, performed for instance, by the OPAL collaboration ||TT 
leads to the constraints on the {Mh, Ma) plane shown in Fig. It follows that the direct 
experimental limits on Mh and Ma are rather weak: the data allow for very light (A^) 
provided A^ (h^) is heavier than ~ 65 (50) GeV even if only tan (3 > 1 is allowed. In particular, 
there exist no absolute bound on Mh from LEP data, provided one respects the bound on 
sin^(/3 — a) shown in Fig. |l]. 

For large values of tan/9, independent constraints on and A^ follow from the non- 
observation of the Yukawa process. The available limits on the {Ma, tan f3) plane are shown in 
Fig. 1^ |jl5[. Similar limits on both (M/j,tan/3) and (M^,tan/?) planes have been reported only 



recently by the DELPHI Collaboration |T6|. 

It should also be mentioned that some, rather weak (and dependent on the assumptions 
made about the Higgs boson decay branching fractions) limits on the very light h° and A° can 
be derived from the so-called Wilczek processes |jl3|, i.e. from T and J/ip decays into h^{A^) 
and the photon [0]. We do not take these limits into account in performing the fit to the 



electroweak data, because the results we will show do not change as M^, varies over the range 
0-30 GeV. It should also be clear that whenever sin^(/9 — a) ~ the heavier CP— even scalar, 
H^, is constrained by the Bjorken process so that Mh ~ 90 GeV. Finally, the non-observation 
of the charged Higgs boson production at LEP sets the bound Mh+ > 72 GeV [O , which is. 



however, much less restrictive than the indirect limit derived in the 2HDM(II) from the 6 — > 57 
process and, for tan /3 < 1, from Rb. 

From the above, it follows that within the 2HDM(II) there still exist two scenarios with 
either very light scalar or very light pseudoscalar A^, which are not excluded by the available 
data (for other constraints not discussed here, see |3^, |33| and Section 5). 



In the next section we will consider how the electroweak precision data constrain these two 
scenarios. 
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Figure 2: 95% CL. limits on the [Mh, Ma) plane from the Bjorken process and the associated 
Higgs boson production as given by the OPAL collaboration |15| . 
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Figure 3: 95% CL. limits on the (tan/3, Ma) plane from the Yukawa process as given by the 
ALEPH collaboration |15|. 
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4 Global fit to the precision data 



Since the advent of LEP and SLAC experiments, precision electroweak data have been playing 
an increasingly important role in constraining the mass of the top quark in the Standard Model 
(SM). Nowadays, with the top quark mass directly measured at Fermilab |jl8[, they significantly 



constrain the last unknown parameter of the SM, the mass of the Higgs boson |19]. They are 
also very useful in constraining possible forms of new physics such as supersymmetry pO] or 



technicolour [^]. It is, therefore, natural to ask to what extent the precision electroweak data 
constrain the Two Higgs Doublet Model of the type II. In this context, a question of particular 
importance is whether they are still compatible with the existence of a very light neutral Higgs 
particle (scalar or pseudoscalar one) which, as we discussed in the preceding Section, is not yet 
excluded by direct searches. In this Section we discuss under what conditions the existence of 
light h'^ or A'^ can be compatible with indirect constraints imposed by precision data. 

A strictly statistical approach to constraining the 2HDM(II) indirectly would consist of 
finding the global minimum of the fit ^ind, in the next step, excluding all points in the 
parameter space for which Ax^ is greater than 3.84 (exclusion at 95% C.L.). Comparison of 
such an analysis with the similar one carried for the SM would reveal that, per degree of freedom 
(d.o.f.), the fit in 2HDM(II) is much worse than in the SM. This follows from the fact that 
the description of the electroweak data by the latter is nearly perfect [|I|. Hence, a 2HDM(II) 
predicting individual observables as accurately as the SM would have much worse x^/d-o.f., 
because of larger number of free parameters in its Higgs potential. In our investigation we do 
not follow such an approach. Rather, we take the SM best global value as a reference point 
and concentrate on the qualitative discussion of which regions of the 2HDM(II) parameter space 
can give an equally good global value. In order not to be too restrictive, when constraining 
the parameter space we use the rough criterion that the in the 2HDM(II) should not be 
greater than the SM best value plus 4. Our emphasis is, however, on the fact that such a 
criterion (in fact any reasonable one) does allow for a very light neutral scalar or pseudoscalar. 
The data we take into account include the precision electroweak data reported at the Moriond 
'98 conference which are dominated by those from LEP 1. For future reference we record 
that the best SM fit to the electroweak data we have chosen gives us ~ 15.5. Therefore, all 
bounds on the masses of the 2HDM(II) Higgs bosons we will present are derived by requiring 
the x^ ill that model to be less than 19.5. 

In discussing the values taken by the x^ for various Higgs boson mass configurations (to 
explain qualitatively the origins of the bounds we show) we will always start with a discussion 
of the contribution to the Ap parameter, defined as 

_ Hw(0) l^zzjO) ^sw ^z^jO) 
^ Ml \^ Ml ' 

where sw [cw) is the sine (cosine) of the Weinberg angle. It largely determines the bulk of 
the predicted values of the electroweak observables such as Mvi/, sin^ d^^^ (measured through 
various asymmetries of the final fermions and/or their polarizations [Q), etc., and is therefore 
the main factor shaping the x^ curves, at least for not too small or too large values of tan/3, 
i.e. when the couplings of the Higgs bosons to the 66 pair are not enhanced. 
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In addition, for low and high values of tan /5, a particularly important role is played in the fit 
by the quantity Rfy. The current experimental result is Rf, = 0.21656, with error ARf, = 0.00074 
(which is 0.9 standard deviation above the SM prediction). Its importance follows from the fact 
that, in the 2HDM(II), the contribution of the Higgs bosons can easily change the prediction for 
Rb with respect to the SM, spoiling the fit to the data. This contribution has been studied 
in detail in ref. |2^. (Handy formulae for 6Rb can also be found in the Appendix of ref. [^.) 
The contribution of H'^ to 6Rb contains parts that are proportional to {rrit / M z)^ cot'^ j3 and 
{irib/Mzy tan^ (3 and, consequently, can be sizeable for either very small or very large values of 
tan/3. The neutral scalars become relevant to Rb only in the latter limit (see. Eqs. (0) and (^)). 
For a qualitative understanding of the results it is sufficient to remember that the contribution 
of is always negative, whereas the contribution of the neutral Higgs bosons can be positive 
provided is not too heavy, say M4 < 100 GeV, and the splitting between its mass and the 
mass of the CP— even scalar or (the one which, for a given angle a, couples more strongly 
to the bb pair) is not too large. For example, in the configuration (typical for the Minimal SUSY 
Standard Model with large tan /?) Mh ~ M4 ^ 70 GeV (and sin^ a 1) the contribution of 
and can easily overcompensate the negative contribution of with mass Mh+ > 100 GeV. 
It is the interplay of the contribution to the Ap parameter and the contribution to 6Rb that is 
responsible for interesting bounds on 2HDM(II) with light (pseudo) scalar particle, which can 
be derived on the basis of the fit- 

Another very important constraint on new physics is the measured value of the branching 
ratio BR{B Xg'y) ||2^. In the context of the 2HDM(II) this measurement can be converted 



into a lower bound on the mass of the charged Higgs boson |2^. Recently a big effort was 



made by various groups to improve the accuracy of the theoretical prediction for this ratio 



2g, Our lower bound on Mjj+, based on 6 — 57 for rrit = 174 GeV, is shown as a 

function of tan/? in Fig. ^ by the solid line. f\ It is (for large tan/5) higher by some 35 GeV 
than the recent estimate [^, which gives Mh± ~ 165 GeV. In view of the well-known exquisite 



sensitivity of the limit on Mf{+ to the details of the analysis, this should be considered as a 
satisfactory agreement. We will, however, try to keep open the possibility that the charged 
Higgs boson can be as light as ~ 165 GeV. 



4.1 Light 

We discuss the light scenario first. For the sake of clarity it is convenient to distinguish two 
mass ranges: i) Mh ^ 30 GeV and ii) Mh > 30 GeV. This division follows from the upper 
bounds imposed by the direct LEP search on the allowed value of sin^(/3 — a) [0 (Fig. |l]): in the 
case i) sin^(/5 — a) < 0.01-0.02 (and for practical purposes can be set to zero); in the case ii) the 

^We compute BR{b sj) with NLO accuracy, following the approach of ref. p9[ s upplemented with 
electromagnetic as well as and 1/™^ corrections setting the parameter S — 0.9 ]28|. The theoretical 

uncertainty is taken into account by computing the rate for fib = 2.4 and 9.6 GeV and then shifting its larger 
(smaller) value upward (downward) by the errors, added in quadrature, related to the uncertainties in as, rrih, 
TTic/rai,, \VtbVf^g/Vcb\'^ , and higher-order electroweak corrections; we do not take into account the variation of the 
scale If the resulting band of theoretical predictions for BR{h 57) has an overlap with the CLEO 95% 
C.L. band, the point is allowed. 



7 



55 
50 
45 
40 
35 
30 
25 
20 
15 



20 40 60 80 

M, (GeV) 



100 



0.216 
0,2155 r 

0,215 
0,2145 Fy' 
of 0.214 - 
0,2135 r 

0,213 
0,2125 b- 
0,212 



20 40 60 80 

M. (GeV) 



100 



Figure 4: and Rb as a function of Mh in the case of triple degeneracy with Md = Mh = 
Ma = Mh+ = 250 GeV and svo?{(3 — a) = 0. Short-dashed, sohd, dot-dashed, long-dashed and 
dotted hues correspond to tan/? =0.5, 1, 5, 20 and 50, respectively. The top mass is fixed to 
= 174 GeV. 

upper bound on sin^(/9 — a) changes roughly linearly (on the logarithmic scale) from ~ 0.02 for 
Mh ~ 30 GeV up to 1 for Mh ~ 90 GeV (see Fig. p. In this paper we will be mainly interested 
in case i) and will essentially not explore the case ii), which requires a more involved analysis. 
It will also prove helpful to consider separately two ranges of the parameter tan/? = f2/fi, 
namely 0.5 < tan/? < 10-20 (small and intermediate) and 20 < tan/5 < 50 (large). Distinct 
properties of these regions follow from the different sensitivity of the important observable 
to the masses of the Higgs bosons for small and intermediate values of tan (3, and for large ones. 

In order to understand the main features of the fit to the electroweak data qualitatively, 
it is instructive to begin the discussion of the light case with the somewhat peculiar limit in 
which the remaining Higgs bosons, iJ^, and , are exactly degenerate and have a common 
mass Mu. In Fig. |^a we show the value of obtained in the 2HDM(II) as a function of Mh for 
different values of tan (3 and Mr, = 250 GeV. Corresponding values of Rb are shown in Fig. ^d. 
In both cases we set sin^(/? — a) = and we keep rrit fixed at 174 GeV. The pattern observed 
in Fig. ^ can be easily understood by checking the contributions to the parameter Ap, which 
can be represented, in the 2HDM [^], in the convenient form: 

= A 2 "^i^A, Mh+) + COs2(/3 - «)Ae + SIT? - «)A„ (8) 

where 

d 

= A 2 .,2 [A{MH+.Mh) - A{Ma, Mh)] + Apsm{Mh) (9) 
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and As is obtained by the exchange ^ Mh- The function 

A{x, y) = A{y, x) = ^x^ + ^y^ - log ^ (10) 

8 8 Ax^ — 

is positive and large ii x ^ y and vanishes for x = y. Finally, 

AP5m(M) = --^^[AiM,Mw)-AiM,Mz)] 



a 



w 



log -VrF - -^T75 ^ log ■ 



(11) 



is the Standard Model Higgs boson contribution to Ap. A good quahty of the fit is obtained 
for 

ApNEw = Ap- ApsMiM,,oJ ~ 0, (12) 

where M^o is a reference SM Higgs boson mass ~ 100 GeV. 

Since in the case of triple degeneracy shown in Fig. ^ Ap is independent of tan/5, sin^(/5 — a) 
and Mh values, for fixed Mo the curves reflect mainly ^ the dependence of Rb on Mh for 
different values of tan/3. For low and moderate values of tan/5 the predicted value of Rb is 
M/j-independent, whereas for large tan/? (> 20), for which the h^bb coupling is enhanced (see 
Eq. (0)), the sensitivity of Rb to Mh becomes crucial. In addition, for light (i.e. small Md, 
~ 200 GeV) its negative contribution to Rb spoils the fit for very small or very large values 
of tan/? (but this effect is M/i-independent). On top of that comes the global sensitivity of 
the fit to Md, which enters through ApsMiMc)- This dependence is exactly the same as the 
dependence of the SM fit (for fixed top quark mass) on the value of the SM Higgs boson mass. 
It follows that the best fit (in the case of triple degeneracy) is obtained for Md as low as possible 
and intermediate values of tan /3. Note that only the case of M^y = 250 GeV and tan /3 = 5 (and 
marginally, for heavier h^, also for tan (3 = 20) shown in Fig. ^ yields an acceptable value of x^- 
Allowing for Md = 165 GeV gives (for intermediate values of tan/3) ^ 17.5. In the region 
i) taking sin^(/5 — a) = 0.01 improves only marginally (it decreases by ~ 0.2 for Mr, = 250 
GeV). In the region ii) taking the largest value of sin^(/3 — a) allowed by the LEP data, e.g. 
sin^(/? -a) = 0.2 for Mh = 50, GeV decreases for Md = 250 GeV the value to ~ 18.1 as a 
result of a "redistribution" of contributions in the last two terms in Eq. (|^). 

Essential improvement of the fit is, however, obtained by the departure from the strict limit 
of triple degeneracy. For very small and intermediate values of tan f3, this is illustrated in Fig. ^ 
where for Mh = 20 GeV (which is representative of < Mh < 30 GeV), sin^(/3 — a) = 0, and 
rrit fixed at 174 GeV, we show x^ as a function of Ma for different combinations of and 
masses. 

To explain the pattern of x^ seen in Figs. |^a for tan/? = 2 (which is representative of 
intermediate values of tan/? for which the predicted value of Rb is as in the SM), recall that 

^For a qualitative explanation of the shapes of the curves the other "oblique" parameters, S and U |]3T| , 
play only a secondary role. 
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Figure 6: a) Lower limits on mass coming separately from b ^ s'y (solid line) and from 
the requirement that Rf^^^^ > R^^^ - 2{ARi,)^^^ (dashed line) as a function of tan/3. Also 
shown are upper limits on Mh+ arising for < 20 GeV from the requirement of fine tuning 
in Ma not larger than 1% and 3% (dotted and dash-dotted hues, respectively), b) Upper hmits 
on Mh in the case of light /i° (M/, <20-30 GeV; sin(^ - a) = 0) and Mh+ = 1000, 800 and 
600 GeV (solid, dashed and dotted lines, respectively). 



the SM-type contribution (pTD to Ap, Eq. (H), is negative and (for sin^(/3 — a) = 0) decreases 
with increasing mass. This too negative contribution can be easily compensated for by the 
contribution of the other Higgs bosons to Ap, when the equality of H^, H'^ and masses is 
relaxed. From Eqs. (||) and (P) it follows that their contribution is positive, provided Mh+ ~ 
Ma- Obviously, the Mh+-Ma mass splitting must increase with increasing Mh in order to 
compensate for the increasingly negative contribution to ApsM- It is also clear that the 
possibility to adjust the total Ap to a proper value should hold also for larger M^ and/or 
sin^(/5 — a) 7^ 0. We have checked that, for example for M^ = 50 GeV and sin^(/3 — a) = 0.2 
(see Fig. |I|), the plots look very similar to those for Mh = 20 GeV and sin^(/5 — a) = 0. 

Since the contribution of H^, and to Ap depends on differences of the masses squared 
(quadratic violation of the 5'f/y(2) "custodial" symmetry), the cancellation between their con- 
tribution and ApsM which is needed in the case of Mh < 20 GeV to adjust Ap to a proper value, 
becomes more and more delicate as the mass of Mh+ increases. This leads to a stronger and 
stronger correlation of Mh+ with Ma, clearly seen in Figs. |^ (note the different mass scales in 
the panels). This correlation becomes particularly fine-tuned in the case of small tan/3, where 
the requirement of good Rf, forces Mh+ to be large. The lower limits imposed on Mh+ by 
b S'-y and Rb are shown in Fig. ^a by the solid and dashed lines, respectively. Since the very 
strong correlation of Mh+ with Ma niay seem unnatural, it is interesting to see how the re- 
quirement of "naturalness" of the constrains the parameter space. This is illustrated by the 
dotted (dash-dotted) line in Fig. ^a, which bounds from below the region in the (tan/3, Mh+) 
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plane in which a change by <1% (3%) of value which gives the minimum of (fo^' that 
particular point in the (tan/3, Mh+) plane), does not lead to > 19.5. In producing these 
curves, the fit was always optimized with respect to the values of and rrit. It is also 
interesting to note (see Figs. ^) that in the case of heavy H'^, only relatively light can givea 
a good fit to the data. Therefore, for fixed values of Mh+ and tan f3 the requirement of 
< 19.5 leads to an upper bound on Mh shown in Fig. (where we have optimized the fit 
with respect to M4 and m^). 

Larger value of at the minimum for heavier observed in Figs. ^ can be explained 



by the behaviour of the parameter S [^, which was not taken into account in the above 
discussion. We define S" at = M| rather than at = (since it is S'(M|) that parametrizes 
more effectively the electroweak observables measured at LEP): 



S 

where 



a 



c'wFzz{Ml) - clF,,{Ml) + ^{2sl - 1)F^,(M|) 

Sw 



(13) 



^ Mil_IhM. ,14) 

For g2 = it is easy to obtain the analytic expression for S{0), which we record in the Appendix. 
Since for approximately fixed the parameter S grows with increasing Mh, as shown in Fig. |^, 
it is obvious that for too high values of Mh the quality of the fit will be spoiled. 

When Mh+ increases, its negative contribution decreases 5*, but this decrease is almost 
entirely compensated by the change in Ma, which is required by the Ap variable so that 
remains in its minimum (for that value of Mh+)- Therefore, the upper limit on Mh is, for 
tan/3 ^ 1.5, almost independent of the assumed value of Mh+ as seen in Fig. ^o. For smaller 
tan j3, however, a lighter induces a larger negative contribution to 6Rb, which has the effect 
that the upper limit on Mh is stronger for lighter H~^, since the minimum of is already 
higher than for heavier H~^. Therefore, this limit is rather stringent for 0.5 < tan/3 < 0.8. 

In the case of light and tan/3 ^ 20, a qualitatively new behaviour of x^ appears as a 
result of the interplay of h^, A'^ and H'^ contributions to (for sin^(/3 — a) ~ the boson 
contributes negligibly to Rb). Because of the tight correlation of Ma with Mh+ (required by 
Ap), no acceptable x^ can be obtained for too heavy since there is then a large negative 
contribution to Rb due to the large mass splitting between and A^. Thus, for large tan (3, very 
light necessarily implies the existence of relatively light H~^. The corresponding pattern of x^ 
(for fixed = 174 GeV) is illustrated in Fig. ^ for two different values of Mh+ and tan (3 = 20, 
35 and 50. Deeper minima of x^ for larger values of Mh seen in Figs, ^c-f can be explained 
by the fact that (for fixed Mh+) larger Mh requires lighter ^4° to give acceptable Ap (the 
correlation seen already in Fig. ^ and this increases the positive contribution to Rb of the 
latter. From the pattern seen in Fig. ||, it follows that there exists, for a given mass of the 

^Whenever we optimize with respect to mt, the top mass measurement m,t = (173. 9±5) GeV is included 
as one of the fitted data in our fit. 
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Figure 7: Parameters 5" as a function of for Mfj+ = 300 GeV and 1 TeV for = (from 
below) 100, 250, 500, 1000 and 2000 GeV. In all cases Mh = 10 GeV and sin2(/? - a) = 0. 

lighter scalar and for a given upper bound on Mh, an upper bound on Mh+. For = 20 
and 10 GeV, this bound (obtained for sin(/? — a) = by scanning over Mh, Ma and nit and 
looking for points with < 19.5) is plotted in Fig. for an assumed upper limit on Mh 
equal to 1000, 3000 and 5000 GeV. Taking /i° lighter strengthens the bound, as can be easily 
inferred from Fig. ^. Of course, for large tan/5 values allowing for heavier weakens the 
upper bound on Mh+- However, for tan/3 values ~20-30, where Rb gradually ceases to be so 
important (even allowing Mh to vary up to 3 TeV), the minimum of (fo^' given Mh+) is 
achieved for Mh = 90 GeV (lower limit on the SM-like Mh from direct searches). Therefore 
the minimum of does not depend on the assumed upper limit on Mh- The importance of 
this upper bound on Mh taken together with the lower one coming from 6 — > 57 is obvious f\. 

Since the best is obtained, for tan P > 25, for the lowest possible mass of the charged 
Higgs boson and Mh equal to its assumed upper bound, for a given value of tan P there exists a 
lower bound on M^ that follows from the assumed upper bound on Mh and the lower limit on 
Mh+ coming from b —>■ sj. This bound (obtained by scanning over rrit and Ma as well as over 
sin^(/5 — a) in the experimentally allowed range shown in Fig. |l|) is shown in Fig. for two 
different lower bounds on Mh+ taken to be 200 GeV (our limit from b —>■ 37) and 250 GeV and 
different upper bounds on Mh- (For Mh+ = 200 GeV, the upper limit on M^ does not decrease 
when larger Mh are allowed, since it would require Ma < 65 GeV, which is excluded by the 

^However, one should note that, for given values of the Higgs sector parameters, the minimum of for 
tan/3 close to 50 occurs for nit ~ 166 GeV, for which the lower bound on Mff+ coming from & ^ 57 is slightly 
(by about 20 GeV) lower than the bound plotted in Fig. |^a (which was obtained for rrif — 174 GeV) 
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Figure 8: as a function of Ma for three different values of tan /3 and different if + masses for 
Mfi — 20 GeV (i.e. sin^(/? — a) = 0). Solid, dashed, dotted, dot-dashed and long-dashed lines 
correspond to Mh = 90, 500, 1000, 2000 and 5000 GeV, respectively. 
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Figure 9: Limits from the fit: a) Upper for the mass as a function of tan/5 for = 20 
(sohd hues) and 10 GeV (dashed hues) assuming the upper hmit on Mh equal (from bottom 
to top) to 1, 3 and 5 TeV. b) Lower for for Mh+ = 200 and 250 GeV (sohd hues) assuming 
the upper hmit on Mjj equal to 1 TeV, and for Mh+ = 250 GeV (dashed line) and Mh < 3 
TeV. 

OPAL analysis, see Fig. |^.) The bound turns out to be very sensitive to the lower limit on 
Mfj+, showing that the b ^ s'j process is crucial for constraining a light scalar Higgs scenario 
for large tan/?. 

The global limits shown in Figs. Pa,b must be confronted with the recent analysis |jl6| of the 
Yukawa process for production. In the case of Fig. ^ parts of the solid (dashed) lines corre- 
sponding to tan /3 > 40 (30) seem to be excluded by the data for the Yukawa process. However, 
those parts of these lines that are not excluded still provide interesting and complementary 
limits on the light scenario. In the case of Fig. the limit obtained for Mh+ = 200 GeV is 
for most of the tan f3 range weaker than the limit imposed by the Yukawa process. For heavier 
these limits shown in Fig. Wp become competitive to the ones derived in ref. [ffH . 



4.2 Light 

In the the case of light the analysis is more involved because there is effectively one more 
variable: sin^(/3 — a). Therefore let us begin again with the case of triple degeneracy, = 
Mh = Mh+- As is clear from Eq. (^ the contribution of the Higgs sector to Ap is in this limit 
independent of the value of Ma and of sin^(/5 — a); it is given simply by ApsM{Mh). Therefore, 
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for M/i = Mh = Mh+ and moderate values of tan/5 (1.5 < tan/3 < 20), the best value of 
depends on the lower limit on Mh+ from b — *• sj. For tan/? > 1 the condition < 19.5 is 
satisfied if Mh+ < 240 GeV, which is still allowed (see Fig. 1). For tan/3 < 1 the light H+ that 
is needed to give good Ap tends to give too negative a contribution to Rb (see Fig. ^) and the 
condition < 19.5 cannot be satisfied. 

Even for Mh+ larger than 240 GeV, the value of can be kept below 19.5 by the departure 
from the limit of triple degeneracy. For example, consider the limit in which and H'^ remain 
degenerate (which makes Ap independent of sin^(/3 — a)). It is then easy to see that for any 
value of Mjj+ (and any Ma) there exists a solution to the equation Ap = that occurs for 
Mh = Mh ~ Mh+- For values of tan/3 < 20 (for which neutral Higgs bosons do not play any 
role in Rb) and Mjj+ not too large (so that the effects of the S parameter are not too large, see 
below) the existence of such a solution is sufficient to ensure small values of x^- Note, however, 
that since the solution to the equation Ap = is due to the cancellation of the Ap^EW, which 
depends on the mass splittings quadratically against ApsM^Mh) which depends on Mh only 
logarithmically, such a solution is strongly fine-tuned (the more, the heavier H~^). The fine- 
tuning can be, however, reduced (or, more precisely, shifted to the variable sin^(/3 — a)) by 
relaxing the condition Mh = Mh- This is illustrated in Fig. |1^ where, for M^ = 10 GeV and 
rrit = 174 GeV, we show as a function of Mh for different (moderate and low) values of tan j3 
and different choices of Mh+ and Mh- In all cases we optimize with respect to the value of 
sin^(/3 — a). 



The peculiar dependence of as a function of Mh seen in Fig. 10 can be understood by 
inspection of the formulae for Ap, Eq. (|^). Consider first the case of Ma = 10 GeV, Mh+ = 300 
GeV, as in Fig. |lO]a. For Mh = 90 GeV and Mh < Mh ^ Mh+, both factors, A^ and A„ 



are positive and obviously cannot cancel the positive first term in Eq. (|^) for any choice of 
sin^(/3 — a). As Mh increases, however, A^ decreases rather fast and takes on negative values for 
Mh ~ Mh+, whilst Ac decreases slowly (logarithmically) and reaches negative values only for 
very large values of Mh- Obviously, for Mh such that A^ < —aA{MA,MH+)/4:nslyM^ < A^ 
there always exists a choice of sin^(/3 — a) for which Ap = 0. This explains the plateau in 
X^ for light and Mh > Mh+- f\ With Mh increasing, Ac decreases very fast. In addition. 
As decreases too, though only logarithmically. As a result, for Mh larger than some critical 
value, both Ac and A^ become negative and smaller than — a74(M^, Mjy+)/47rs^M^, and 
again Ap cannot vanish leading to large values of x^ seen in Figs. p!0|a-c for mass configurations 
corresponding to dot-dashed and long-dashed lines. It should also be obvious that for very 
heavy the fine-tuning in sin^(/3 — a) becomes extremely big, making such solutions rather 
unnatural. Thus light A'^ and values of tan/3 < 1, for which Mh+ must be large, are rather 
unlikely. 

For very large Mh+, Ap can vanish only for heavy H^- In the case of heavy this means 
that the parameter S is also large (see Fig. |^) and spoils the quality of the fit. This explains 
why, in Fig. |TO|c, for larger values of Mh, and even for the other masses in configurations for 
which Ap can vanish by a judicious adjustment of sin^(/3 — a), the value of x^ is still above 

^Eventually, for Mh large enough, also Ac becomes smaller than —aA{MA,MH+)/4:^Ts'^,M^,, so that there 
is again no solution with Ap = 0. This happens for smaller Mh, the smaller the value of Mh+- 
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Figure 10: ^-s a function of Af// for different low and intermediate values of tan/3 and 
different masses, for = 10 GeV. Solid, dashed, dotted, dot-dashed and long-dashed 
lines correspond to equal (90, 150, 200, 250, 275) GeV, (90, 200, 400, 450, 475) GeV and 
(90, 250, 500, 750, 940) GeV for panels a, b and c, respectively. 
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19.5. 



From the above explanation it is clear that for small and intermediate values of tan (3 a 
light CP— odd scalar can be tolerated provided is lighter than some bound Mb, which is 
only slightly smaller than the mass of the charged Higgs boson (which in turn is constrained by 
6 — *• S7 and, for tan/5 < 1, by the Rf, measurements). At the same time, Mh is bounded from 
below also by roughly the same mass Mb- Of course, according to our previous discussion, the 
case Mh = Mb < Mh+ is always allowed (in this case sin^(/5 — a) is completely unconstrained). 
It is also worth noting that, for M^, <C Mb f» Mb ~ Mb+, the value of sin^(/3 — a) (which is 
crucial for the h'^ and production processes) that is needed to keep below 19.5 is close 
to 1, implying that Mh is constrained, in this case, by the LEP search to be greater than ^ 90 
GeV. For heavier H^, sin^(/? — a) decreases, which means that the experimental lower bound 
on Mh is also relaxed appropriately (recall that the limits from the associated h^A^ production 
require only Mh Z 50 (70) GeV for M4 = 10 (50) GeV - see Fig. |^ - irrespectively of the value 
of sin^(/? — «)). The decrease of sin^(/3 — a) with Mb > Mb+ is slightly faster for heavier AP 
(and/or h^) and slower for heavier . 

For large tan (3 the dependence of for different values of Mh and Mb+ is shown 



in Fig. |TT|. The mass of the CP— odd scalar is taken to be 25 GeV so as to respect the bound 
from the Yukawa process also for tan/5 = 50 [|1^, |T6|. The behaviour of can be explained by 
combining the information contained in Figs. (reflecting mainly the behaviour of Ap) with 
the behaviour of the corrections to Rh for different Higgs boson mass configurations. In order 
for corrections to the latter quantity not to be too negative the neutral scalar that couples more 
strongly to the bb pair cannot be too heavy. On the other hand, for Mb ~ Mb+, a good Ap is 
obtained for sin^(/3 — a) ~ 1 (i.e. sin a ?s 0), which means that it is the mass splitting between 
and A^ that is relevant to Pf,. Only for sufficiently heavy P° does sin^(/3 — a) become small 
enough (sin a 1) so that couples with full strength to 66, and becomes relevant to P^. 
Thus, a good fit to the data can be obtained only with light and rather heavy P° (since 
this occurs for sin^(/5 — a) <^ 1, can be lighter than 90 GeV as we have just explained). 
The increase of for Mb > 1 TeV seen in Fig. |ll] a and d for Mh = 70 and 90 GeV is due to 
the fact that, for lighter P+, Ac becomes smaller than —aA{MB+ , Ma) /4:^TS^^M^r already for 
relatively light H^. Since, as explained above, Ac decreases very fast as Mh increases, this effect 
is even more pronounced for Mh = 150 GeV. Note that because of this effect for Mb+ = 200 
GeV one can reach < 19.5 only for Mh < 90 GeV. 

For the same value of Mb and Mh the is slightly larger for heavier (despite the 
fact that the negative contribution of the latter Higgs boson to Rf, is decreased), because of 
the behaviour of sin^(/3 — a): it is larger for heavier and therefore, light does not fully 
compensate for the effects of light ^4° as couples more weakly to the bb pair. The conclusion 
following from the above considerations and from Figs. |Tl| is that, in the scenario with large 
tan/3 and light CP— odd neutral Higgs particle {Ma ~ 25 GeV), the mass of the lighter neutral 
CP— even boson h'^ is bounded from above. This upper bounds on Mh for four different values 
of the light CP— odd scalar mass (obtained by scanning over rrit, sin(/? — a), Mb and Mb+) 
are shown as functions of tan /? in Fig. Of course, as tan (3 decreases, the upper limit on 
Mh approaches the bound Mb ~ Mb+ discussed previously. The existence of an upper bound 
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Figure 11: as a function of Mh for tan/? = 35 and 50 and different i/+ masses. Solid, 
dashed, dotted lines correspond to Mh equal to 70, 90 and 150, respectively. Ma = 25 GeV. 
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Figure 12: Upper limit on in the light CP— odd scalar scenario as a function of tan/? for 
Ma = 25 GeV (sohd line), 15 GeV (dashed line), 10 GeV (dotted line) and 1 GeV (dot-dashed 
line) . 

on Mfi means that LEP2 will be able to effectively test the scenario with light A° and large 
tan/? either via the Bjorken process or the associated Higgs boson production, at least for 
tan/? > 45(30) for Ma < 25 (10) GeV. 



5 Summary of the results of the global fit and other 
constraints 

We have investigated the impact of the precision electroweak data on the parameter space of 
the 2HDM of type II. We have been particularly interested in constraints imposed on the very 
light neutral scalar scenarios (with hP or ^4° in the range < 25-30 GeV) by the requirement of 
good (as good as in the SM) fit to the electroweak data. It turns out that neither scenario is 
excluded (directly or indirectly), provided some constraints are respected. 

Apart from the well-known constraints on the charged Higgs boson mass set in 2HDM(II) 

by 6 ^ S7 and, for tan/? < 1, by Rb, we find that in the case of a light CP— even scalar 
and of intermediate values of tan/? masses of the CP— odd and charged Higgs bosons must 
be tightly correlated (the more, the heavier is H~^) in order to maintain in the 2HDM(II) the 
same quahty of the fit to the data as in the SM; this leads to strong "fine-tuning" of this 
scenario. Therefore, limiting the acceptable degree of "fine-tuning" yields an upper bound on 
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Mu+ (Fig. I^a). For a given mass of the requirement of a good fit to the data puts also an 
upper hmit on Mh-, which is particularly strong for tan/5 < 1 (Fig. ^b). For large values of 
tan P (tan P > 20-30) the interplay between the corrections to Ap and implies that if hP is 
light, must also be light. The upper limit (as a function of tan/?) depends on the assumed 
upper bound on Mh but, for Mh in the TeV range, it is very strong for tan/? > 35, and Mh+ 
of the order of 300 GeV (Fig. ||a). It is also interesting that, for large tan/3 and Mh+ > 200 
GeV, the electroweak data set the lower limit on the mass of shown in Fig. ^d. 

For the light scenario we find that, for low and intermediate values of tan/?, the quality 
of the fit is maintained provided M^ < Mb < Mh, where Mb ^ Mh+- For tan/? > 30-35 there 
emerges an additional upper bound on the mass of the lighter scalar which enables the test 
of the light A°-large tan/3 configuration at LEP2. 

Other potential sources of further constraints on the 2HDM(II) that we have not considered 
here are the muon g — 2 measurement and Z — > h{A)'y decays p5|] . 

The upper limit on BR{Z^ — > h^{A^)'-y) decays set by LEPl constrain our model only for 
h^{A^) masses below 20 GeV and for values of tan/3 either very low (below 0.2) or very large, 
above 55-70 GeV (which we have not considered). For larger (A°) masses, tan/? is pushed 
outside the range (|^) in which perturbative calculations can be done reliably. In this case, the 
resulting constraints can eventually be given a meaning as constraining the effective (on-shell) 
Z^h^'^ and Z^A^'^ couplings but the relation of these couplings to the original parameters of 
the model (which we have been using) cannot be calculated perturbatively. 

The constraints following from the present measurement of the g—2 of the muon are stronger 
than the ones following from the Yukawa process only for /i°(y4°) masses below 1-2 GeV, in 
which case they exclude values of tan/3 > 4 (10) for Mh(A) = 0.1 (1) GeV. The E821 experiment 
may (depending on the measured central value) improve these limits (according to the analysis 
presented in [0), so that they become stronger than the Yukawa process ones up to M^^a) ~ 10 
GeV and exclude values of tan/? > 2 (15) for M^^a) = 1 (10) GeV. 



Appendix 



Here we give the expression for the parameter S'(O). It differs slightly from the formula presented 



in ref. 36 



S = -^ log ^ + - sm\f3 - a) \A'{Ma, Mh) + A'{Mw, M,) + MIB'{Mz, M,) 

OTT Mw TT L 

+ - cos2(/3 - a) \A'{Ma, Mh) + A'{Mw, Mh) + MIB'{Mz, Mh 



The functions A' and B' are given by: 
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